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Abstract 



We present data on multiple jet production for transverse energies greater 
than 20 GeV in pp collisions at y/s = 1.8 TeV. QCD calculations in the parton 
shower approximation (pythia) and in the next-to-leading order approxima- 
tion (jetrad) show discrepancies with data for three and four-jet production. 
This disagreement is especially apparent in multiple jet angular and transverse 
momentum distributions. 
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The study of jet production at high transverse energy was one of the main goals of 
the 1993-1995 run of the Fermilab Tevatron colhder, and the results have been compared 
with leading-order QCD predictions by both the CDF |I| and D0 collaborations. These 
high-ii^T data, where i?T is the transverse energy of the jet were described satisfactorily by 
complete tree-level leading order 2 —>■ N QCD calculations and by the herwig parton 
shower Monte Carlo [|]. In this paper, we describe studies of the complementary kinematic 
region of Q'^/s <^ 1, where is the square of the momentum transfer between partons, 
which we set equal to E^, and s is the square of center of mass energy in the rest frame of 
the collision. Here the BFKL ||^ description of jet production differs significantly from that 
of the high-ii^x DGLAP ||^ kinematic domain of ~ s. Measurement of jet production in 
this kinematic region can provide information on the evolution of higher-order jet processes. 

We present results that extend our previous measurements of multiple jet production to 
lower Et^. The data were collected with the D0 detector during 1993-1995 at a proton- 
antiproton center-of-mass energy of 1800 GeV. Jets were measured in the liquid-argon 
calorimeter, which has a segmentation of Arj x A0 = 0.1 x 0.1, where t] is pseudorapidity 
and (p is azimuthal angle At least one calorimeter trigger tower [At] x A0 = 0.2 x 0.2) 
with Et > 2 GeV was required by the Level-1 trigger, and at least one jet with E^ > 12 
GeV was required by the Level-2 trigger Jets were reconstructed using a fixed cone 
algorithm with radius ATI = y/Ar]'^ + = 0.7 in 77 — space. The jet reconstruction 
threshold was 8 GeV. If two jets overlapped and the shared transverse energy was more 
than 50% of the transverse energy of the lower-energy jet, the jets were merged; otherwise 
they were split into two jets. The integrated luminosity of this data sample was 1.96±0.29 
nb~^. Instantaneous luminosity was restricted to be below 3 xlO^°cm~^s~^ to minimize 
multiple pp interactions. 

To provide events of high quality, online and offline selection criteria were used to suppress 
multiple interactions, cosmic ray backgrounds, and spurious jets. Jets were restricted to the 
pseudorapidity interval \rj\ < 3. 

Jet energies have been corrected for calorimeter response, shower development, different 
sources of noise, and contributions from the underlying event [^. These corrections comprise 
the largest source of systematic uncertainty on the jet cross section. The typical value of the 
correction to jet energy is 15 - 30%, with an uncertainty of 2-4%. In our study, we consider 
jets with E't > 20 GeV. For an n-inclusive jet event, the n leading jets must have transverse 
energy above the threshold value. The trigger efflciency is 0.85 for the inclusive (n = 1) jet 
sample for energies near threshold, rising rapidly to unity at larger i?T- The efficiency is 
essentially unity for n > 1. 

To compare with data, Monte Carlo (MC) events were generated using the pythia 6.127 
||T0| and JETRAD [0] programs. These generators simulate particle-level jets in the parton- 
shower approximation, and parton-level jets in the next-to-leading order approximation, for 
PYTHIA and JETRAD, respectively. The smearing of jet transverse energies was implemented 
using the experimentally determined jet energy resolution 0, which is ~ 20% at E^^ = 20 
GeV. In PYTHIA, jets were reconstructed at the particle level using the D0 algorithm, and 
in JETRAD, at the parton level, using the Snowmass algorithm |Q. 

Distributions in transverse energy for the leading jet for n=l to n=4 inclusive jet events 
are shown in Fig. |l], together with the results from PYTHIA simulations. In these and all 
other plots, the data has been corrected for inefficiencies and energy calibration, but not 
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for contributions from an underlying event. Also, we normalize the theory (increased by a 
factor of 1.3) to the observed two-jet cross section Fig. |l](b) for i?T > 40 GeV. Figure |^ 
shows the fractional difference (Data - Theory) / Theory for the i?T spectra in Fig. |I] with the 
systematic uncertainties, arrising from uncertainties in jet-energy calibration and resolution. 
The theory is in agreement with the data for the single-inclusive jet sample in the entire 
interval, and with the two-jet sample for most of the energy interval (some excess of data 
is observed at low Et). However, for the three and four-jet samples, there is large excess 
relative to theory at low i?T and a deficit near 75 GeV. The shapes of the experimental and 
theoretical spectra are clearly different, and not reconcilable through re- normalization. 

The systematic uncertainty on the cross section is due primarily to the uncertainty in the 
energy calibration. The uncertainty from energy resolution represents the main uncertainty 
in the MC. The uncertainty from the energy calibration can be estimated by considering 
spectra with ±1 standard-deviation corrections to jet E^^. The same procedure can be 
used to derive the uncertainty due to jet resolution in the MC. At 25 GeV uncertainty in 
the three-jet cross section due to calibration is 36%, and the uncertainty in the MC due 
to resolution is 17%. In Fig. ^ the relative systematic uncertainties corresponding to the 
energy calibration added to a 15% uncertainty in luminosity are shown (in quadrature) by 
the solid lines centered about zero. The uncertainties from energy smearing are shown by 
the dashed lines near the data points. The total systematic uncertainties on the ratio are 
shown by the dotted lines. Because the systematic uncertainties are highly correlated in 
Et (a change of the cross section in one bin is accompanied by a corresponding change 
in neighboring bins), the departure of the ratio from zero cannot be explained solely by 
systematic uncertainties. 

To explore the discrepancies in three and four-jet production, we turn to observations of 
azimuthal distributions, distributions in summed transverse momenta, and three-jet studies. 
In Fig. |^(a) we plot the azimuthal difference between the leading two jets in events with 
two or more jets. Figures ^(b)-|^(d) show the azimuthal difference between the first and sec- 
ond, first and third, and second and third highest-i?T jets in a three-jet event. In Fig. |^(a) 
we see the strong anticorrelation (in the transverse plane) expected of two-jet events. The 
distribution widens substantially in the three-jet sample (Fig. |^(b)-|](d)). The peaks corre- 
spond to the kinematic constraint of transverse momentum conservation for jets produced 
in hard QCD subprocesses. Altough, in general, pythia reproduces the observed shapes, 
there is a large excess of events in the three-jet sample not consistent with expectation. In 
particular, there is a significant contribution to three-jet events with two jets back-to-back 
in the transverse plane near ($ = vr). 

Distributions in the square of the summed vector transverse momenta of jets = 
(Eti + Et2 + ■ ■ ■ + Etw)^ shown in Fig. H(a-c) indicate that the excess corresponds to events 
with a large imbalance. In fact, when these events are removed (with a requirement of a good 
balance in transverse momentum), the three-jet data and theory come into better agreement 
at small E^. The shoulder at Q"^ ~ 1600 GeV^ in Fig. ^(a) is eliminated when the event 
sample is restricted to just two jets with E^ above 20 GeV, and no other jets between 8 and 
20 GeV. This shoulder can consequently be associated with higher-order radiation. 
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FIG. 1. Distributions in the transverse energy of the leading jet for (a) single-inclusive, (b) 
two-jet inclusive, (c) three-jet inclusive and (d) four-jet inclusive events. Histograms show the 
PYTHIA simulation normalized (increased by a factor of 1.3) to the inclusive two-jet sample at 
> 40 GeV. 



To find the pair of jets {i,j} most likely to originate from the hard interaction (rather 
then from gluon Brehmsshtrahlung) , we define the scaled summed dijet vector transverse 
momenta: q^j = (Eti + ETj)/{ETi + Etj). We choose the pair with the smallest magnitude 
of this vector and plot the distribution of the relative azimuthal angle $c between the jets 
in that pair Fig. ^(a). The data lie above theory in the region where two jets, reflecting 
a hard scatter, appears back-to-back ($c = tt). PYTHIA shows a broader distribution, and 
the prediction from JETRAD is peaked away from $c = ti" due to the presence of the third 
(radiated) jet. 
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FIG. 2. (Data - Theory)/Theory as a function of the transverse energy of the leading jet for 
(a) single-jet inclusive, (b) two-jet inclusive, (c) three-jet inclusive and (d) four-jet inclusive event 
samples. 



Figures |](b) and Fig. ^(c) show the azimuthal separation of the third jet from each of 
the two jets that correspond to the minimum qfj. These distributions contain events only 
for TT — $c < 0.4; that is, events in which the balanced jets are essentially back-to-back. 
When the third jet is correlated with the balanced jets, it will be expected to be emitted 
along or opposite to the balanced jets. The uncertainties from the energy calibration and 
luminosity are shown by the solid lines, and from the energy resolution by dashed lines. We 
see that the data has a wider distribution than pythia, and much wider distribution than 
JETRAD. The third jet appeares to be uncorrelated with the balanced jets, and is emitted 
at all angles. Our studies indicate that the observed differences in shape are not sensitive 
to modeling of the underlying event or contributions from multiple-parton scattering. 
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<I>j3 (rad) (> 3 jets) (^ad) (> 3 jets) 

FIG. 3. Distributions of the relative azimuthal angle between two jets in (a) two-jet events 
and in three-jet events (b-d). Jets are ordered by their transverse energies. Histograms are from a 
PYTHIA simulation of such events. 

In summary, our data on multiple-jet production at low show significant discrepancies 
with PYTHIA and jetrad. This is observed in the distributions of the transverse energy of 
the leading jets (Fig. P, in the square of the summed vector transverse momenta Q'^ (Fig. §), 
and in the three-jet angular distributions that suggest the presence of an uncorrelated jet 
(Fig. Additional corrections to QCD calculations are therefore required to accommodate 
these results; higher-order or BFKL processes are possible candidates. 

We thank the staffs at Fermilab and collaborating institutions for their contributions 
to this work, and acknowledge support from the Department of Energy and National Sci- 
ence Foundation (U.S.A.), Commissariat a L'Energie Atomique (France), State Committee 
for Science and Technology and Ministry for Atomic Energy (Russia), CAPES and CNPq 
(Brazil), Departments of Atomic Energy and Science and Education (India), Colciencias 
(Colombia), CONACyT (Mexico), Ministry of Education and KOSEF (Korea), and CON- 
ICET and UBACyT (Argentina). 



9 




2000 4000 2000 4000 2000 4000 



(GeV') (> 2 jet) (GeV') (> 3 jet) Q^2 (QeV") (> 4 jet) 

FIG. 4. Distributions in the square of the summed vector transverse momenta Qj^, for two-jet 
inclusive, three-jet inclusive and four-jet inclusive event samples (a-c). Histograms show the pythia 
simulation. 




(rad) (> 3 jets) (rad) (> 3 jets) (^ad) 3 jets) 

FIG. 5. The distribution of the relative azimuthal angle in three-jet events, between the jets 
in the pair with the minimal scaled summed transverse momentum (a), and between the third jet 
and the other two leading Et jets in the pair (b-c) . Histograms show the PYTHIA simulation, and 
the open symbols the JETRAD simulation. 
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